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1. Introduction 

The exclusive production of mesons was studied in detail mostly close to 
the kinematical threshold. The Tevatron is a first accelerator which opens a 
possibility to study the central (semi) exclusive production of mesons at high 
energies. A similar program will be carried out in the future at the LHC. 
Here I review several mechanisms of exclusive meson production studied 
recently by our group (the details can be found in [U [2J, [3J, SI O El E] ) - 
In general, the mechanism of the reaction depends on the quantum num- 
bers of the meson and/or its internal structure. For heavy scalar mesons 
(scalar quarkonia, scalar glueballs) the mechanism of the production, shown 
in FigfTJ is exactly the same as for the diffractive Higgs boson production 
extensively discussed in recent years by the Durham group [9]. The dom- 
inant mechanism for the exclusive heavy vector meson production is quite 
different. Here there are two dominant processes shown in Figj2j When 
going to lower energies the mechanism of the meson production becoming 
more complicated and usually there exist more mechanisms. For example in 
Figj3]I show a new mechanism of the glueball production proposed recently 
in Ref.0. 

Other exclusive or semi-exclusive processes were discussed during the 
conference by A. Martin [TO], Ch. Royon [11] and W. Guryn |12j . 



* Presented at the Epihany 2009 conference: "Hadron interactions at the dawn of the 
LHC", 5-7 January 2009; dedicated to the memory of Jan Kwiecihski. 
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Fig. 1. A sketch of the bare QCD mechanism of exclusive heavy scalar meson 
production. 




Fig. 2. Two basic QED ® QCD mechanisms of exclusive heavy vector meson pro- 
duction. 

2. Selected examples 

Recently we have calculated differential cross sections for several exclu- 
sive processes: 



• pp - 


-* ppr]' , pp — * ppr] c ( IP IP + 77 


• pp - 


-» pp X c(0 + ) (IP IP + 77 ) 


• pp - 


- PPX c(l + ) (IP IP) 


• pp - 


-> pp/ (1500) (IP IP + 7T+7T-) 


• pp - 


-> ppJ/tfj (IP 7 + 7 IP) 
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Fig. 3. A sketch of the bare QCD mechanism of exclusive heavy scalar /o(1500) 
meson production. 



• pp — > ppT ( IP 7 + 7 IP) 

• pp -» ppir+ir- ( (IP + IR) ® (IP + IR)) 

• AA —* AAp°p° (77) 

Above the dominant mechanisms are shown in the parantheses. 

The details of the formalism as well as a detailed analysis of differential 
distribution in longitudinal and transverse momenta can be found in our 
original papers [TJ [21 El El El E] ) . Here I wish to discuss only some illustrative 
examples. 



2.1. Exclusive Xc production in proton-proton and proton- antiproton 

collisions 

According to the Khoze-Martin-Ryskin approach (KMR) [9], we write 
the amplitude of the exclusive double diffractive color singlet production 

pp -> ppxcj as 

fg{ f (xu x[, ql p q\ p t 1 )fg{ f (x 2 ,x' 2 ,ql t , q\ t M) 

— 2 — T^i • V 1 

%,t 1\,t lit 

The amplitude is averaged over the color indices and over the two transverse 
polarizations of the incoming gluons [9J . 

In calculating the vertex Vj lC2 we have included off-shellness of gluons 
[3]. The unintegrated gluon distributions were taken from the literature. We 
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have demonstrated in Ref.[3] that for relatively light Xc(0), unlike for the 
Higgs boson [9], the dominant contributions come from the nonperturbative 
regions of rather small gluon transverse momenta. 

In Ref.[3] we have made a detailed presentation of differential distri- 
butions. Here only selected results will be shown. As an example I show 
distribution in Feynman variable xf of the Xc meson for three different en- 
ergies: W = 200 GeV (RHIC), W = 1960 GeV (Tevatron) and W = 14000 
GeV (LHC) for different UGDFs from the literature. Characteristic for 
central diffr active production all distributions peak at xf ~ 0. Although 
all UGDFs give a similar quality description of the low- a; HERA data for 
the i<2 structure function, they give quite different longitudinal momentum 
distributions of Xc(0 + ). The UGDFs which take into account saturation 
effects (GBW, KL) give much lower cross section than the BFKL UGDF 
(dash-dotted line). Therefore the process considered here would help, at 
least in principle, to constrain rather poorly known UGDFs. 




Fig. 4. Distribution in Feynman x F for RHIC (W=200 GeV), Tevatron (W=1960 
GeV) and LHC (W=14000 GeV) for different UGDFs: BFKL (dash-dotted), KL 
(dashed), GBW (dotted) and KMR (solid). 

The three-body reactions lead to correlations of outgoing protons. In 
Figj5]I show an example for diffractive mechanism with KL UGDF (as an 
example) as well as for photon-photon fusion [3J. The {t\,t%) distribution 
obtained in the photon-photon fusion mechanism differs qualitatively from 
the distribution of the diffractive mechanism. One can see a strong en- 
hancement of the cross section when ti, — ► which is caused by the 
photon propagators. Although the diffractive component is subjected to 
much stronger absorption effects than the electromagnetic one, it is clear 
that the diffractive component dominates. 

The situation with the axial-vector production is new compared to both 
zero-spin case (scalar [3], pseudoscalar [1] mesons) as well as to the vec- 
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Fig. 5. Two-dimensional maps in t\ and t% for KL UGDF (left) and for two-photon 
fusion (right). 

tor meson production where the vector meson is dominantly transversely 
polarized [21 [5], at least for small transferred four-momenta in the nucle- 
onic line. The axial-vector meson can be polarized both transversely and 
longitudinally. 

There is interesting theoretical aspect of the double diffractive produc- 
tion of the Xc(l + ) meson. The coupling ggx c {^- ++ ) vanishes for on-shell 
gluons (so-called Landau- Yang theorem) . According to the original Landau- 
Yang theorem [13] the symmetries under space rotation and inversion forbid 
the decay of the spin-1 particle into two (on-shell) spin-1 particles (two pho- 
tons, two gluons). The same is true for the fusion of two on-shell gluons. 
The symmetry arguments cannot be strictly applied for off-shell gluons. 

In Ref . [1] we have confirmed explicitly that the Landau- Yang theorem is 
violated by virtual effects in diffractive production of Xc(l + ) leading to very 
important observational consequences. In our approach the off-shell effects 
are treated explicitly. For comparison, in the standard KMR approach the 
corresponding cross section would vanish due to their on-shell approxima- 
tion. The measurement of the cross section can be therefore a good test of 
the off-shell effects and consequently UGDFs used in the calculation. 

In Fig. 01 show distributions in rapidity y for different UGDFs from the 
literature. The results for different UGDFs differ significantly. The biggest 
cross section is obtained with BFKL UGDF and the smallest cross section 
with GBW UGDF. The big spread of the results is due to quite different 
distributions of UGDFs in gluon transverse momenta (qu,Q2t), although 
when integrated over transverse momenta distributions in longitudinal mo- 
mentum fraction (2:1,0:2) are fairly similar. 

Comparing the left and right panels, the cross section for the axial- 
vector x c (l + ) production is much smaller (about two orders of magnitude) 
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Fig. 6. Distribution in rapidity of Xc(l + ) meson (left panel) and x c (0 + ) meson 
(right panel) for different UGDFs. 



than the cross section for the scalar Xc(0 + ) production. This is related to 
the Landau- Yang theorem, which "causes" vanishing of the cross section 
for on-shell gluons. For axial-vector quarkonia the effect is purely of off- 
shell nature and is due to the interplay of the off-shell matrix element and 
off-diagonal UGDFs. This interplay causes a huge sensitivity to UGDFs 
observed in Fig. [6j 

At the Tevatron the \c mesons are measured through the 7 + J/ty decay 
channel. The axial- vector Xc(l + ) meson has a large branching fraction for 
radiative decay Xc(l + ) — ► 7 + J/ip (BR = 0.36 [H]). This is much bigger 
than for the scalar Xc(0 + ) where it is only about 1 % [J3]. Therefore the 
discussed off-shell efects are very important to understand the situation in 
the 7 + J/ty channel observed experimentally. 



2.2. Exclusive /o(1500) production 
in proton-proton and proton- antiproton collisions 

In Ref.fB] we have discussed exclusive production of scalar /q(1500) in 
the following reactions: 

P + p^P + Ml500)+p, (2) 
P + p^P + Ml500)+p, (3) 
p + p^ n + / (1500) +n . (4) 



While the first process could be measured at the J-PARC complex be- 
ing completed recently, the latter two reactions could be measured by the 
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Fig. 7. The integrated cross section as a function of the center of mass energy for 
pp — > pp/o(1500) (left panel) and pp — > nn/ (1500) (right panel) reactions. The 
thick solid lines are for pion-pion MEC contribution (A = 0.8, 1.2 GeV), the dashed 
line is for QCD diffractive contribution obtained with the Kharzeev-Levin UGDF, 
the dotted line for the KMR approach and the thin solid lines (blue on-line) are 
for "mixed" UGDF (KL ® Gaussian) with a = 0.5, 1 GeV. The dash-dotted line 
represents the two-gluon impact factor result. 



PANDA Collaboration at the new complex FAIR planned in GSI Darm- 
stadt. The combination of these processes could shed more light on the 
mechanism of /o(1500) production as well as on its nature. 

In Ref.[6j we have proposed a new mechanism (see Figj3|) which becomes 
dominant at lower energies. In FigJ7]we show the integrated cross section for 
the exclusive /o(1500) elastic production pp — > p/o(1500)p and for double 
charge exchange reaction pp — > n/o(1500)n. The thick solid line represents 
the pion-pion component calculated with monopole vertex form factors with 
A = 0.8 GeV (lower) and A = 1.2 GeV (upper). The difference between 
the lower and upper curves represents uncertainties on the pion-pion com- 
ponent. The pion-pion contribution grows quickly from the threshold, takes 
maximum at W ~ 6-7 GeV and then slowly drops with increasing energy. 
The gluonic contribution calculated with unintegrated gluon distributions 
drops with decreasing energy towards the kinematical threshold and seems 
to be about order of magnitude smaller than the pion-pion component at W 
= 10 GeV. We show the result with Kharzeev-Levin UGDF (dashed line) 
which includes gluon saturation effects relevant for small-x, Khoze-Martin- 
Ryskin UGDF (dotted line) used for the exclusive production of the Higgs 
boson and the result with the "mixed prescription" (KL <g> Gaussian) [6] for 
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different values of the do parameter: 0.5 GeV (upper thin solid line), 1.0 
GeV (lower thin solid line) . In the latter case results rather strongly depend 
on the value of the smearing parameter. 



2. 3. Exclusive production of T 
in proton-proton and proton- antiproton collisions 

The photoproduction amplitude is the major building block for our pre- 
diction of exclusive T production in hadronic collisions. The amplitude for 
the reaction under consideration is shown schematically in FigEJ As it is 




Fig. 8. A sketch of the exclusive jp — > Tp amplitude. 



explained in Ref . [15] , the imaginary part of the amplitude for the 7* p — > Tp 
process can be written as 



9m M x „ Xv (W,t = -A 2 ,Q 2 ) = W ^^^ J *£ as tf)F{x uX2 , KuK2 ) 

/dzd^k 
r x _ n I\~,,\ v (z, k, Ki, K 2 , Q 2 ) ,(5) 

where the transverse momenta of gluons coupled to the QQ pair can be 
written as 

A A 

«i = « + y , k 2 = -K + — . (6) 

The quantity J 7 (x\,x 2 , k%, k 2 ) is the off diagonal unintegrated gluon distri- 
bution. Explicit expressions for I\ 1: \ v can be found in [15]. For heavy vec- 
tor mesons, helicity-flip transitions may be neglected, and we concentrate 
on the s-channel helicity conserving amplitude, A 7 = Ay. In the forward 
scattering limit, i.e. for A = 0, azimuthal integrations can be performed 
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analytically, and we obtain the following representation for the imaginary 
part of the amplitude for forward photoproduction jp — > Tp : 



Qm M {W, A 2 = 0, Q 2 = 0) = W 2 c ^^ a ™ 



4vr 2 



1 dz 



z(l - z) 



irdk 2 ipy(z, k 2 



K 



-as 



(q 2 )f(x eff , k 2 )(a (z, k 2 ) W (k 2 , k 2 ) + k 2 ) W x (k 2 , k 2 ) 



where 



and 



A (z,k 2 ) 
Ax{z,k 2 ) 

W (k 2 ,K 2 ) 
W x (k 2 ,K 2 ) 



2 , k 2 m b 
mi + 



M + 2m b 
z 2 + (1 - z) 2 - (2z - l) 2 



M + 2m 6 



k' 2 



k 2 + m 



2 ■ 



0) 
(10) 



k 2 + m 2 / ^2 - m 2 - k 2 ) 2 + 4ra 2 fc 2 



A; 2 + m 2 
2k 2 



1 + 



k 2 — m 2 — k 2 



(k 2 — m 2 — k 2 ) 2 + Am 2 k 2 
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We treat the T, T' mesons as 66 s-wave states, the relevant formalism of 
light-cone wave functions is reviewed in |15j . 

The necessary formalism for the calculation of amplitudes and cross- 
sections was outlined in detail in Ref. [2]. Here I give only a brief sum- 
mary. The basic mechanisms are shown in FigJUJ The major difference from 




Fig. 9. A sketch of the two mechanisms considered: photon-pomeron (left) and 
pomeron-photon (right), including absorptive corrections. 
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photoproduction at ep machines, where the photon was emitted by a lepton 
which does not participate in the strong interactions, now, both initial state 
hadrons can be the source of the photon. Therefore, it is now necessary 
to take account of the interference between two amplitudes. The photon 
exchange parts of the amplitude, involve only very small, predominantly 
transverse momentum transfers. Here we concentrate on the kinematic do- 
main, where the outgoing protons lose only tiny fractions z±, z 2 <C 1 of their 
longitudinal momenta, in practice z < 0.1 means y < 3. In terms of the 
transverse momenta of outgoing hadrons, p 1 2 , the relevant four-momentum 
transfers are U = — (pf + zfrn^)/{l — z%) ,i = 1, 2, and si « (1 — Z2)s and 
s 2 ~ (1 — z\)s are the familiar Mandelstam variables for the appropri- 
ate subsystems. Photon virtualities Qf are small (what counts here is that 
Q\ <C My), so that the contribution from longitudinal photons can be safely 
neglected. Also, as mentioned above, we assume the s-channel-helicity con- 
servation in the 7* — > T transition. 

The 2^3 Born-amplitude (without absorptive corrections) can be writ- 
ten in the form of a two-dimensional vector (corresponding to the two trans- 
verse (linear) polarizations of the final state vector meson): 

M(°)( Pl ,p 2 ) = e 1 -^F KXl (p 1 ,t 1 )M rh2 ^ V h 2 (s2,t2,Q 2 1 ) + (1 <- 2)2) 

Z\ tl 

Inclusion of absorptive corrections (the 'elastic rescattering') leads in mo- 
mentum space to the full, absorbed amplitude 

M( Pl ,p 2 ) = J ^S el (k)M^(p 1 -k,p 2 + k) = M^(p 1 ,p 2 )-SM(p 



With 

S el (k) = (27r)V 2 ) (k) - \T{k) , T(k) = aZ(s) exp ( - \B el k 2 ) , (14) 

where o 3 ^ t {s) = 76 m b, B e \ = 17 GeV~ 2 were taken, the absorptive correc- 
tion 5M reads 

6M( Pl ,p 2 ) = J J^T(k)MM( Pl -k,p 2 + k). (15) 
The differential cross section is given in terms of M as 

da = * 2 |M| 2 dydhdhdcf) , (16) 
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Fig. 10. o~tot{lP — * T(lS')p) as a function of the jp cm-energy versus HERA-data. 
Left: dependence on the treatment of the bb — > T transition; solid curves: Gaussian 
(G) wave function, dashed curves: Coulomb-like (C) wave function. Thick lines 
were obtained including the NLO-correction for the T decay width, while for the 
thin lines Knlo = 1- Right: dependence on the slope parameter Bq (given in 
GeV -2 ), for the Gaussian wave function. The experimental data are taken from 

puna ni 



where y is the rapidity of the vector meson, and <p is the angle between p 1 
and p 2 . 

In FigllQII show the total cross section for the exclusive "fp — ► Tp process 
as a function of the -fp cm-energy. In the left panel I show results for two 
different wave functions discussed in the text: Gaussian (solid lines) and 
Coulomb-like (dashed lines). Free parameters of the wave function have 
been adjusted to reproduce the leptonic decay width in two ways: (a) using 
leading order formula (thin lines) and (b) inlcuding QCD corrections (thick 
lines). Including the -ff 7VLO _iac t or i n the width enhances the momentum- 
space integral over the wave function (the WF at the spatial origin), and 
hence enhances the prediction for the photoproduction cross section. The 
ratio of the cross section for the first radial excitation T(2S) to the cross 
section for the ground state T(1S) is shown in Figllli The principal reason 
behind the suppression of the 2S state is the well-known node effect - a 
cancellation of strength in the 2S case due to the change of sign of the 
radial wave function. It is not surprising, that the numerical value of the 
2S'/l,S'-ratio is strongly sensitive to the shape of the radial light-cone wave 
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Fig. 11. The 2S/lS'-ratio er tot ( 7 p -> T(2S)p)/a tot (jp -> T(1S» as a function of 
the 7p cm-energy. 



In our calculations we assumed an equality of the slopes for T(15) and 
T(2S) production. This appears to be justified, given the large spread of 
predictions from different wave functions. We finally note, that the ratio 
depends very little on the choice of the K^lo factor (compare left and right 
panel). 



2.4- Exclusive production of the tt + tt pairs 
in proton-proton collisions 

Up to now I have studied only exclusive production of a single mesons. 
Also the channels with meson pairs seem interesting. In particular, the 
channel with two charged pions which seem feasible experimentally. 

The underlying mechanism was proposed long ago in Ref.[20j. The gen- 
eral situation is sketched in Fig|12i The corresponding amplitude for the 
pp — ► ppir + ir~ process (with four-momenta p a +Pb — ► Pi +P2+P3 +P&) can 
be written as 

M*™™ = M 13 (h,s 13 ) F{t a ) — F(t a ) M 24 (t 2 ,s 24 ) 

t a ~ mi 

+ M 14 (ti, 8 U ) F(t b ) F{t b ) M 23 (t 2 , 813) , (17) 
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Fig. 12. A sketch of the dominant mechanisms of exclusive production of the ir + ir 
pairs at high energies. 

where denotes "interaction" between nucleon i=l (forward nucleon) or 
i=2 (backward nucleon) and one of the two pions k = tt + (3), tt~ (4). In 
the Regge phenomenology they can be written as: 




Above sik = Wf k , where Wik is the center-of-mass energy in the (i,k) subsys- 
tem. The first terms describe the subleading reggeon exchanges while the 
second terms describe exchange of the leading (pomeron) trajectory. The 
strength parameters of the ttN interaction are taken from Ref.|21j. More 
details of the calculation will be presented elsewhere [7]. The 2 — > 4 am- 
plitude (fT7|) is used to calculate the corresponding cross section including 
limitations of the four-body phase-space. 

Here I wish to show only one example of the two-dimensional distribution 
in rapidity of positively charged pion and rapidity of negatively charged 
pion at the LHC energy of W = 14 TeV. The distribution observed differs 
considerably from the distribution of the phase space factor. One can see a 
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Fig. 13. Rapidity distribution of tt + versus it for W = 14 TeV. 



two-dimensional shape of the ridge form elongated along the line 2/3 = 2/4- 
The minimum of the cross section on the top of the ridge occurs when 2/3 = 
y4 = and two maxima close to the phase space ends. The minimum occurs 
in the part of the phase space where the pomeron-pomeron contribution 
dominates, i.e. when both Wik are comparable and large. The maxima 
are related to the dominance of the pomeron-reggeon and reggeon-pomeron 
mechanisms, i.e. where one of is small and the second one is large. 
The reggeon-reggeon contribution is completely negligible which is due to 
the fact that both Wik cannot be small simultaneously. We hope that the 
ALICE collaboration at the LHC will be able to measure such distributions. 



2.5. Exclusive AA -> AAp°p° 
in ultrarelativistic collisions 

Exclusive production of elementary particles (lepton pairs, Higgs, etc.) 
or mesons (vector mesons, pair of pseudoscalar mesons, etc.) in ultrarel- 
ativistic collisions is an interesting and quickly growing field |22| [23l [24"] 
of theoretical investigation. On experimental side the situation is slightly 
different. So far only single-/? exclusive cross section AA — > AAp° was mea- 
sured [25]. Here the dominant mechanism is a photon-pomeron (pomeron- 
photon) fusion. 

Let us consider the process AA — ► AAp°p° depicted in Fig]14[ The cross 
section takes the familiar form of a convolution of equivalent photon fluxes 
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Fig. 14. The reaction discussed in this paper. 



and 77-cross sections: 

da(AA -> p°p°AA; saa) , , , ~ , 00 \ 
= dn 77 (xi,x 2 ,o)<T(77 ->■ p p ;x 1 x 2 saa) + • -( 19 ) 

The effective photon flux is expressed through the electric field strengths of 
the ions [8]. 

Often flux factors of equivalent, almost on-shell, photons are calculated 
as for point-like particles with rescaled charge e — ► Ze, and the total cross 
section is calulated using a simple parton-model type formula: 



a (AA -> A(p°p )A) = J du^^^a (77 - P°P°) • (20) 



n(a>i) n(w 2 ) 

The formulae (|20p clearly does not take into account absorption effects when 
initial nuclei undergo nuclear breakup. This can be easily done in the impact 
parameter space where the geometry of the collision is more explicit. Then 
rather two-dimensional flux factors [26] must be used. 
The simple EPA formula can be generalized to 

a (AA -» A(p°p )A) = J d 2 b 1 daj 1 d 2 b2dLo 2 N(uj 1 ,b 1 )N(L02,b2) 

6 (|6i - b 2 \ - R12) a (77 -> p°p°) . (21) 

Here an extra 6 function was introduced which excludes those cases when 
nucler collisions, leading to nuclear breakup, take place (R12 = Ri + -^2)- 
The two-dimensional fluxes in (|2ip can be calculated in terms of the charge 
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form factor of nucleus |28j as: 

N(w,b) = ^*{x,b), (22) 
where the auxiliary function $ reads: 

•(«,»)- r*,^(.)ft^j3ffl'. (2 3 ) 

The second ingredient of our approach is 77 — ► p°/9° cross section. Here 
the situation is not well established. The cross section for this process was 
measured up to W 77 = 4 GeV [29]. At low energy one observes a huge 
increase of the cross section. 

In Fig|15l we have collected the world data (see [29] and references 
therein). We use rather directly experimental data in order to evaluate 
the cross section in nucleus-nucleus collisions. In Fig(15] we show our fit to 
the world data . 



* TASSO 
□ CELLO 

O TPC/2-GAMMA- 

♦ PLUTO 
ARGUS 



— VDM-Regge : 




Fig. 15. The elementary cross section for the 77 — > p°p° reaction. In the left panel 
we display the collection of the e + e~ experimental data [5S] and our fit. In the 
right panel we show our predictions based on the VDM-Regge model decribed in 
the text. For comparison we show also result when the form factor correcting for 
off-shell effect is ignored (see [8]). 

The cross section above W = 4 GeV was never measured in the past. It 
is well known that the cross section for 77 — ► hadrons can be well described 
in the VDM-Regge type model. We use a similar approach for the final state 
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channel p°p°. In Figll5lwe present the corresponding i-integrated cross sec- 
tion together with existing experimental data taken from [29] . The vanishing 
of the VDM-Regge cross section at W 77 = 2m p is due to t m i m t max limita- 
tions. It is obvious from Fig {15] that the VDM-Regge model cannot explain 
the huge close-to-threshold enhancement. In Fig(16]we show distribution of 
the cross section for the nucleus-nucleus scattering in photon-photon center- 
of-mass energy for both low-energy component and high-energy VDM-Regge 
component. Below W = 2 GeV the low-energy component dominates. The 
situation reverses above W = 2 GeV. One can study the high energy com- 
ponent imposing extra cut on M pp . However, the cross section drops quickly 
with increasing invariant mass of two-/) mesons. 




Fig. 16. The Au + Au — > Au + Au + p°p° cross section as a function of W 77 = M pp 
for the RHIC energy y/s NN = 200 GeV. The low- and high-energy components are 
shown separately. 

For illustration in Fig(T7] we show the model distribution in impact pa- 
rameter b = \b% — &2 1 . We show distributions for the low- and high-energy 
components separately. I also show distributions for point-like charge, 
monopole form factor and realistic charge density (see [8]). One can see 
slightly different results for different approaches how to calculate flux fac- 
tors of equivalent photons. 

Finally in Fig[18]I show distribution in rapidity of the p°p° pair. Com- 
pared to the point-like case, the distribution obtained with realistic charge 
density is concentrated at midrapidities, and configurations when both p°'s 
are in very forward or both p 0, s are in very backward directions are strongly 
damped compared to the case with point-like nucleus charges. One can also 
see a difference between results obtained with an approximate monopole 
form factor and with the exact one calculated from realistic charge density. 
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Fig. 17. The Au + Au — ► Au + Au + p°p° cross section as a function of the impact 
parameter b for \/s NN = 200 GeV. The meaning of the curves is the same as in 
FigHU The cut off for R12 « 14 fm is clearly visible. 




Fig. 18. The Au + Au — ► Au + Au + p a p a cross section as a function of the rapidity 
of the p°p° pair Y for \/~s NN = 200 GeV. The meaning of the curves is the same 
as in FigfPol 



3. Conclusions 

We have derived, for the first time in the literature, QCD amplitude 
for exclusive elastic double diffractive production of axial-vector x c (l + ) me- 
son. According to the Landau- Yang theorem the amplitude vanishes for 
the fusion of on-shell gluons. We have generalized the formalism proposed 
recently for diffractive production of the Higgs boson and derived corre- 
sponding g*g* — » Xc(l + ) vertex function. The predicted total cross section, 
obtained from the bare amplitude is of the order of a fraction of nb, de- 
pending on the model of UGDF. This is about two orders of magnitude less 
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than a similar cross section for Xc(0 + ). 

However, because the branching fraction BR{x c {^ + ) — ► J/ip + 7) 3> 
Si?(Xc(0 + ) — > J/i/' + t), one may expect a different situation in the J/V> + 7 
channel. This is similar as for inclusive production of P-wave quarkonia, 
where the signal (in the J/ip + 7 channel) of Xc(l + ) is of similar size as that 
for Xc (0 + ). 

We have calculated the forward amplitude for 7p — > Tp reaction within 
the formalism of /^-factorization. In this approach the energy dependence 
of the process is encoded in the x-dependence of unintegrated gluon distri- 
butions. The latter object is constrained by data on inclusive deep inelastic 
scattering. The i-dependence for the jp — > Tp process involves a free pa- 
rameter and is parametrized. We have used different Ansatze for the bb 
wave functions. The results for T(15) production depend only slightly on 
the model of the wave function, while the 2S/1S ratio shows a substantial 
sensitivity. We compared our results for the total cross section with a re- 
cent data from HERA. Our results are somewhat lower than data, although 
the overall discrepancy is not worrysome, given the large uncertainties due 
to the rather poor experimental resolution in the meson mass. The ampli- 
tudes for the 7J» — ► Tp process are used next to calculate the amplitude for 
the pp — > ppT reaction assuming the photon-Pomeron (Pomeron-photon) 
underlying dynamics. We have calculated several differential distributions 
including soft absorption effects not included so far in the literature. Our 
predictions are relevant for current experiments at the Tevatron. 

For the first time in the literature we have estimated the cross section 
for exclusive /o(1500) meson (glueball candidate) production not far from 
the threshold. We have included both gluon induced diffractive mechanism 
and the pion-pion exchange contributions. 

The first component was obtained by extrapolating down the cross sec- 
tion in the Khoze-Martin-Ryskin approach with unintegrated gluon distribu- 
tions from the literature as well as using two-gluon impact factor approach. 
A rather large uncertainties are associated with the diffractive component. 
At present only upper limit can be obtained for the diffractive component 
as the /o(1500) — > gg decay coupling constant remains unknown. The cou- 
pling constant could be extracted only in high-energy exclusive production 
of /o(1500) where other mechanisms are negligible. 

The calculation of the meson-exchange contribution requires introducing 
extra vertex form factors. At largest PANDA energies they are relatively 
well known and the pion-pion fusion can be reliably calculated. We predict 
the dominance of the pion-pion contribution close to the threshold. Our 
calculation shows that the diffractive component is by more than order 
of magnitude smaller than the pion-pion fusion component in the energy 
region of the future PANDA experiments. The diffractive component may 
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dominate over the pion-pion component only for center-of-mass energies 
W > 15 GeV. Clearly an experimental program is required to disentagle 
the reaction mechanism. 

We have made first estimate of the exclusive production of pairs of 
charged pions at high-energy. Different combinations of pomeron-reggeon 
fusion were included. Rather large cross sections are predicted at RHIC, 
Tevatron and LHC. Here I have shown only rapidity distributions of pi- 
ons at LHC. The production of the two pions is strongly correlated in 
the (y(7r + ), y(ir~)) space. Even at relatively high energies the inclusion 
of reggeon exchanges is crucial as amplitudes with different combination of 
exchanges interfere or/and irN subsystem energies can be relatively small 
W n N < 10 GeV. At high-energies we find a preference for the same hemi- 
sphere (same-sign rapidity) emission of tt + and tt~ . For example at LHC 
energies the same hemisphere emission constitutes about 90 % of all cases. 

We have calculated, for the first time, realistic cross sections for exclusive 
p°p° production in ultrarelativistic heavy- ion collisions at RHIC. We have 
used realistic charge densities to calculate the nuclear charge form factors. 
The absorption effects have been included. 

It was shown that calculating both flux factors and 77 — ► p°p° realisti- 
cally is necessary to make reliable etimates of the nucleus-nucleus exclusive 
production of the p°p° pairs. Large cross sections, of the order of fraction 
of milibarn, have been found. The bulk of the cross section is, however, 
concentrated in low photon-photon energies (low p°p° invariant masses) 

The p° mesons, decaying into 7r + 7r~ , can be measured e.g. by the STAR 
detector at RHIC and the ALICE detector at LHC. 
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